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Abstract 



McVittie's solution of Einstein's field equations, representing a point mass 
embedded into an isotropic universe, possesses a scalar curvature singularity 
at proper radius R = 2m. The singularity is space- like and precedes, in the 
expanding case, all other events in the space-time. It is shown here that this 
singularity is gravitationally weak, and the possible structure of the region 
R < 2m is investigated. A characterization of this solution which does not 
involve asymptotics is given. 
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I. INTRODUCTION 



McVittie's solution of Einstein's field equations M arises uniquely under the following 
hypotheses 0: 

(i) Space-time (M, g) is spherically symmetric and is filled with a shear-free perfect fluid. 

(ii) g ~ g as R — > oo, where i? is the proper radius function of (M,g) and g is the line 
element of a Robertson- Walker (RW) universe. 

(iii) The energy density of (M,g) is spatially homogeneous; p = pit). 

In two previous papers, we studied various properties of space-times described by (i)- 
(iii). In P| (paper I), we proved the uniqueness result referred to above and discussed basic 
properties of the resulting space-times. (The plural here refers to the two different cases 
which arise according as the RW background has flat (k = 0) or negatively curved (k = — 1) 
spatial sections. Since (ii) above refers to infinite spatial distances, the case k = +1 was not 
dealt with. That gap will be filled here.) In || (paper II), we studied the global structure 
of McVittie's k = solution [jl], which is the unique space-time arising in the case k = 0. 
The line element of this space-time may be written as 

ds 2 = -( l ~ m l 2u \ 2 d f + e m)(i + m / 2u )\dr 2 + rW), (1.1) 
\ 1 + m/2u J 

where u = re 13 ^ 2 , m is constant and duj 2 is the line element of the unit 2-sphere. Notice that 
the case m = gives a k = RW universe, while f3 =constant gives the isotropic form of the 
Schwarzschild solution. The coordinate ranges must be carefully considered, for the radius 
function 

R(u) = u(l + m/2u) 2 (1.2) 

is double valued on the interval u G (0, oo). In fact we showed in paper II that u G (0, m/2] 
corresponds to R G [2m, oo), yielding a solution which represents a point mass embedded 
in a k = RW universe with scale factor S 2 (t) = e _/3( -'\ while u G [m/2, oo) corresponds to 
the same with scale factor S 2 (t) = . We deal with the latter case, and using the proper 
radius R as a coordinate, the line element ( |1 . 1| ) may be written 

rfs 2 = _ (1 _ 2m _ l p 2 R 2 )dt 2 _ _ ^-l/2 dRdt + (1 _ ^)-l rfjR 2 + R 2 duJ 2 ^ 3) 

R 4 R R 

The subscript indicates partial differentiation with respect to t. We refer to the space-time 
with this line element as McV^. 

Notice then that the line element is singular at R = 2m. This is not a coordinate effect, 
as the form of the energy density and pressure (calculated via Einstein's field equations) 
show: 

Sttp = \dl 8n P = - 3 -fi - /3«(1 - ^)" 1 / 2 . (1.4) 
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Thus there is a scalar curvature singularity at R = 2m. The expansion of the fluid flow lines 
is homogeneous, and equal to that of the RW background; 

0(t) = \m- 

Thus we may refer to a given McV^ as expanding or contracting according as 9 is positive 
or negative. In paper II, we showed that in very general circumstances, namely those 
which lead to a big bang (big crunch) for the expanding (contracting) RW background, the 
surface R = 2m acts as a space-like barrier which precedes (succeeds) all other events in the 
space-time. More precisely, we showed that all radial null geodesies originate (terminate) 
at R = 2m. Thus this singular surface shuts off the region R < 2m from the rest of the 
space-time. 

In fact, the coordinates of neither (|1.1| ) nor ( |1.3|) include the region R < 2m. Consider 



the analogy of the Schwarzschild solution. Using isotropic coordinates, one can only access 
the region R > 2m of the space-time. But then the coordinate transformation (|1.2|) puts 
the line element into a form in which the extension into the region R < 2m arises naturally. 
This involves skipping over a coordinate singularity at R = 2m, but nonetheless produces 
the line element for the region R < 2m. Here, the corresponding line element ( |1.3| ) does 
not apply to the region R < 2m. We must ask therefore if the solution does exist in this 
region, and if so, how the space-time metric and curvature extend through the singularity 
at R = 2m. 

Central to these questions is the issue of the gravitational strength of the singularity 
There does not seem to be any way to extend through a sigularity which is classified 
as gravitationally strong (see below for what is meant by this); however we will show that 
the singularity at R = 2m is gravitationally weak, and thus the possibility of extending the 
space-time through this surface remains open. 

We then ask what form the space-time region R < 2m might take. We provide a 
partial and negative answer to this question by showing that the hypotheses (i)-(iii) above 
may be cast in a form which avoids any global assumptions (apart from that of spherical 
symmetry), i.e. (ii) above is replaced with a local statement. It then arises immediately 
that these equivalent conditions cannot be satisfied in the region R < 2m. In this sense, 
McVittie's solution cannot be extended into the region R < 2m. This local characterization 
may be used to describe a point mass embedded in an arbitrary (open, flat or closed) RW 
universe. 

We show that if the local characterization is applied to a tachyon fluid, then a solution 
does exist in, and only in, the region R < 2m. 

We conclude with some comments on the significance of these results. 



II. GRAVITATIONAL STRENGTH OF R = 2M 

Two ways of characterizing the strength of a singularity are in common usage: the 
limiting focusing condition of Krolak || and the strong curvature condition of Tipler [|]]. We 
will use the latter. This is based on the idea of classifying a singularity as strong if it crushes 
to zero volume any extended body which falls into it. This is modeled in the following way 
@. Let 7 : [r , 0) —>■ M be a time-like geodesic which is incomplete at the parameter (proper 
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time along 7) value r = 0. Let Z%*\ a = 1, 2, 3 be linearly independent Jacobi fields defined 
along 7 which (a) are orthogonal to k a , the unit tangent to 7 and (b) vanish at r = T\ > r . 
These define a volume element along 7, given by V(r) = \Z^> A Z^ 2 ) A -Z^f. Then 7 is 
said to terminate at a strong curvature singularity if lim r __,o ^( r ) = for all possible choices 
of Z^ and T\. A similar characterization exists for null geodesies involving two space-like 
Jacobi fields. Then a space-time singularity p is called gravitationally strong if every causal 
geodesic terminating at p terminates at a strong curvature singularity. The singularity 
is called gravitationally weak if the volume elements all have finite non-zero limits at the 
singularity. 

We show that every radial causal geodesic terminating at R = 2m terminates at a weak 
curvature singularity, and so every point on R = 2m is gravitationally weak. We apply some 
general results on singularities in spherical symmetry 0. 

Consider the time-like case. By linearity of the geodesic deviation equation, which is 
satisfied by every Jacobi field, and the definition of V(r) which involves the norm of a 3- 
form, V(t) may be invariantly described in terms of the norm a of a radial space-like Jacobi 
field and the norms X^ A \ A = 1, 2 of space-like Jacobi fields tangent to the metric 2-spheres 
0. The result is 

V(x) = aX^X®. 

For a radial null geodesic, we obtain the same form for V with a set equal to one. 

It can be shown that the Jacobi conditions for X^ A \ for both the null and time-like cases, 
lead to 



XM(t) = R(r) f 

J T 



dv 



n R 2 (vy 

where X^ a \ti) = (this is part of the definition of the Jacobi field). Thus in the null case 
and for a non-central (R > 0) singularity, the weak curvature condition is always satisfied, 
and in the time-like case, V — > if and only if a(r) — > 0. The Jacobi conditions for a include 

a + F(r)a = 0, (2.1) 

where an overdot indicates covariant differentiation along the geodesic in question and 

F{x) = 4* 2 + 2-§ - ™ 

Here, ^2 is the Newman-Penrose Weyl tensor term calculated on a null tetrad based on 
the principal null directions of the spherically symmetric space-time (in which sense it is 
invariant), E is the Misner-Sharp gravitational energy and Ric is the Ricci scalar. See 
for details. The behaviour of a in the limit r — > is thus determined by the limiting 
behaviour of F along the geodesic approaching R = 2m. It only remains to determine this 
behaviour for an arbitrary time-like geodesic. 
For MCV4, we have 

F( r ) - _ _ A-1/2 

1 > ~ i?3 4 2 ' 
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where A = 1 — 2m/ R. Thus the factor determining the limiting behaviour of a(r) as 
R — > 2m + along the geodesic is A(r). By definition, lim T ^ — 0. What is important is 
the rate at which this approach to zero occurs. This is determined by the geodesic equations. 
It turns out that this information can be obtained quite easily. 
The first integral of the radial time-like geodesic equations is 

-(A - ^{3 2 R 2 )i 2 - fall A" 1 / 2 iR + A^i? 2 = -I. 

Solving for t, we can write 

-(A - \ftR 2 )i = ~/3 t i? A' 1 / 2 R ± (R 2 + A - \(3 2 R 2 ) l/2 . 

This equation must always have a solution, and so the first integral reveals the constraint 
(substituting for R in terms of A) 

4m 2 (l - A)- 4 A 2 + 4 A -m 2 (3 2 (l - A)' 2 > 0. (2.2) 

We emphasize that this is not a supplementary condition which we must impose on the 
functions i?(r),t(r), but rather one which necessarily arises as a consequence of the other 
geodesic equations. A similar situation arises in Schwarzschild's space-time; there the con- 
straint is 

m 2 (l- A)" 4 A 2 + A > 0. 
In the limit r — > 0, we know that A — > 0, so that 

lim A 2 > lim -8 2 

Along any particular geodesic, A(t) remains finite in the approach to the singularity. This 
is a consequence of a result of paper II mentioned above, that R = 2m cuts off the big 
bang (crunch) (3 t — > oo from the rest of the space-time. Also, (3 t ^ 0, as we assume a 
non-static situation. Thus either A = 0(1) as r — > 0, or A becomes infinite in the limit, 
but with a finite integral. Therefore we must have A = 0(t p ) for some p, < p < 1. Then 
A -1 / 2 = 0(r~ p / 2 ), and the range of p then indicates that 

lim t 2 F(t) = 0. 

This last result is sufficient to show that, quite generally, a(r) is finite and non-zero in 
the limit r — > 0, yielding a weak curvature singularity along the geodesic in question. See 
for details, and || for the general theory. 

Thus all Jacobi fields have finite norm as the singularity is approached along any radial 
causal geodesic, and so the singular surface R = 2m is classified as being gravitationally 
weak. 
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III. A LOCAL CHARACTERIZATION OF THE SOLUTION 



We have shown that the singular surface R = 2m is gravitationally weak. In paper 
II, we showed that it is space-like and trapped, and so using a triad of Jacobi fields at 
each point along the surface, it may be possible to use this as an initial data surface from 
which to evolve Einstein's equations into the region R < 2m. We defer investigation of 
this possibility and instead attempt to solve Einstien's equations directly, subject to a set 
of conditions equivalent to (i)-(iii) above, but which apply throughout space-time. We find 
that it is impossible to do so. 

Condition (i) above leads to the line element (see section 14.2 of M) 



ds 2 = -e^dt 2 + e^ l \dr 2 + r 2 du 2 ), 



subject to 



/ - \f - V = F(r)e^ 2 , (3.1) 



where the prime indicates partial differentiation with respect to r. Spatial homogeneity 
(condition (iii)) is equivalent to F(r) = 6m/r 3 for some constant m, which turns out to be 
the Schwarzschild mass of the vacuum limit. This is also equivalent to the renormalized 



Hawking mass |1(| of the space-time being equal to constant m. See paper I for details. 



This form of F allows us to obtain the first integral of the second of (|3T 

rR> = R(l-^ + A(t)R 2 ) 1 / 2 , (3.2) 
R 

where R = re M//2 is the proper radius. Then condition (ii) is satisfied if and only if the 
arbitrary function of integration A(t) is equal to zero. 

Leaving condition (ii) aside for the moment, we note that the function A(t) appears in 
the energy density and the Misner- Sharp energy of the space-time: we have 

8np(t) = ^-3A(t), 
4 

E = -irR 3 p(t) + m. 
3 

So neither of these terms isolates A(t). However, it turns out that 

g ab D a RD b R = 1 - — + A(t)R 2 , 
ri 

where D a = (g a b + u a Ub)V h is the spatial gradient orthogonal to the fluid flow, which has 
unit time-like tangent u a = e~ u ^ 2 5^. This quantity, and the terms m and R on the right 
hand side are covariantly defined, and so we have a covariant definition of A(t). So we can 
replace the asymptotic condition (ii) with the local condition 

V 777 

(ii)' g ab D a RD b R = 1 - (3.3) 
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Since the spatial gradient of any quantity is space-like, we see that we must have R > 2m. 
Thus we have demonstrated the following result: 

Proposition: There is no equivalent of McVittie's solution in the region R < 2m. 

The local characterization above applies to the case k = 0, but may be generalized to k 
positive and negative. In the case k = — 1, we showed in paper I that subject to conditions 
(i) and (iii) above, the appropriate value for the function of integration A(t) generated by 
an asymptotic condition is e - ^*-*, where S(t) = e^l 2 is the scale factor of a RW universe 
with negatively curved spatial sections. This leads to 

g ab D a RD b R = 1 - ^ + e^R 2 . 

On the other hand, in an arbitrary RW universe we have 

g ab D a RD b R = 1 - ke^R 2 . 

Thus we can change (ii)' so that it applies to the embedding of a point mass in any RW 
universe: 

(ii)" g ab D a RD b R = 1 - t!IL - ke^R 2 . (3.4) 

R 

Then conditions (i), (ii)" and (iii) describe the space-time representing a point mass m 
embedded in a RW universe with curvature index k and scale factor e^l 2 . 

Unfortunately, the solution in the cases k ^ cannot be obtained in closed form, as ( |3.4|) 
leads to an elliptic integral. Nevertheless, we were able to demonstrate in paper I several 
of the properties of the k — — 1 solution. Among these was that the vacuum limit yields, 
as one would expect, the Schwarzschild solution. A serious doubt is cast on the point mass 
interpretation for the k = +1 solution by the fact that this limit does not exist in this case. 
The energy density of the space-time is given by 

87rp(t) = |/3 t 2 + Zke~^ 

which for k = +1 cannot be set equal to zero, as it can for k = and k = —1, in which 
cases the pressure also automatically vanishes, yielding the spherical vacuum with non-zero 
Weyl tensor, i.e. Schwarzschild space-time. The inability to obtain the Schwarzschild limit 
is equivalent to the fact that flat space-time is not a special case of a k = +1 RW universe. 
Indeed this fact indicates that one cannot embed the Schwarzschild solution into a k = +1 
RW universe in any physical way. In such a situation, one must be able to identify in a 
gauge invariant manner the Schwarzschild part of the solution, which would involve looking 
for a well defined vacuum limit. If one could produce the Schwarzschild solution in this way, 
then by setting m = one would also be able to produce flat space-time, which, as we have 
just seen, is impossible. 

IV. A TACHYON FLUID IN R < 2M 

Two facts suggest that a tachyon fluid might be an appropriate matter distribution for 
the region R < 2m. 
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Firstly, what hinders the existence of the solution in R < 2m is the fact that one of the 
conditions we are trying to satisfy is g ab D a RDi,R — 1 — where D a is the spatial gradient 
orthogonal to the fluid flow. Obviously if this were a temporal gradient orthogonal to a unit 
space-like vector, the situation would be reversed; the equation could only be satisfied in 
the region R < 2m. 

Secondly, in the Schwarzschild solution, the preferred vector field d/dt is time-like in the 
region R > 2m, but space-like in R < 2m. This suggests that something similar may occur 
in McV^, the preferred vector field - the eigen- vector of the energy- momentum tensor whose 
eigen-value has multiplicity one - switches over from being time-like to being space-like. 
We investigate this possibility here, and show that the space-time which arises uniquely is 
remarkably similar to McV^. 

Thus we consider a space-time described by conditions (i)-(iii) below. 

(i) Space-time (M, g) is spherically symmetric and filled with a shear-free tachyon fluid. 

The condition on the energy-momentum tensor is thus 

Tab = (p - P)n a n h + pg ab , 

where g a bn a n b = 1, and we refer to p and p as the tachyon energy density and pressure. The 
shear of a congruence of space-like curves (the tachyon flow lines) is defined in exactly the 
same way as that of a congruence of time-like curves. We can use co-moving coordinates 
and write the line-element as 

ds 2 = e u dt 2 + e^-dr 2 + r 2 du 2 ), 

where v = u(r,t), /i = p{r,t) and the shear- free condition has been incorporated. The unit 
tangent to the tachyon flow is n a = e~ u ^ 2 5^. The vanishing of non-diagonal terms in T a b and 
the tachyon pressure isotropy lead to 



^ = ^1 (4.1) 



1 ,9. 1 



H" - -p' 2 - V - - = F(r)e^ 2 . (4.2) 
2 r r z 

(ii) The tachyon energy density is homogeneous; D a p : = (g a b ~ n a n b)V b P — 0. 
This is equivalent to p = pit), and as in the fluid case, leads to the equation F = —c/r 3 



for some constant c which will be identified later. Then (|4.2| ) may be integrated once to 
obtain 

g ab D a RD b R = 1 - t^cR- 1 - B(t)R 2 . 

The term B(t) appears in the tachyon energy density; 

8np = ^e~ 9(t) - 3B(t). 

By analogy with the fluid case, we take B = to be the third condition. This condition 
guarantees that the c = background (c will be identified as being proportional to the 
Schwarzschild mass) is an isotropic tachyon fluid with Lorentzian sections t =constant of 
zero curvature. The condition could be modified to include constant non-zero curvature, 
but we wish to focus on the case analogous to k — 0. 
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(iii) The gradient of the radius function orthogonal to the tachyon flow obeys g ab D a RDbR = 
1-fciT 1 . 

This last equation may be integrated to obtain 

e M/2 = iL(i _|_ S i n ( a (A _|_ ]nr)), 
3r 

where a(t) is arbitrary. We can then calculate 

e" = 4a t 2 (*)e^(!£-l). 

Rescaling t allows us to set Aa 2 (t)e 9 ^ = 1, and we define a(t) = 2a (t). Carrying out the 
transformation R(r, t) = re M//2 puts the line element into the form 

ds 2 = (-1 + |C _ 1 2 R 2 )dt 2 + ati J(|| _ l)-^dRdt " (|| " l)" 1 ^ 2 + 

3i? 4 3i? 3i? 

The vacuum limit must give the Schwarzschild solution. We can show that this limit is 
a =constant, which gives the line element 

ds 2 = (-1 + |^ 2 - - + R *du\ 

and so we identify c = 3m, where m is the Schwarzschild mass parameter. Then the line 
element reads 

= ( _ x + 2m _ 1 2 R 2 )dt 2 + atjR (^ _ i)-V*dRdt - - l)' 1 ^ 2 + (4.3) 
The tachyon energy density and pressure are given by 

3 q 3 o , 2777- . i /o 

87rp = -a t , 8vrp = -a t + a tt {— - 1) 7 . 

The solution only exists in the region R < 2m. The limit m = yields the unique 
isotropic tachyon fluid solution with Lorentzian sections orthogonal to the tachyon flow of 
zero curvature. The vacuum limit a = constant yields the Schwarzschild solution. Setting 
p + p = in McVi, i- e - considering the special case of an Einstein space, we obtain the 
Schwarzschild-de Sitter universe. It is interesting to note that for the tachyon fluid, we 
obtain the Schwarzschild-anti de Sitter universe. For T a b oc g a b is equivalent to p — p = 
for the tachyon fluid, corresponding here to at = A =constant. The metric tensor may then 
be diagonalized using a transformation T = t + f(R), with / chosen to remove off-diagonal 
terms. The resulting line element is 

ds * = _ (1 _ 2m + ^ R 2 )dT 2 + (1 _ 2m + #_ R 2yi dR 2 + ^2^2 
R 4 R 4 

This is the Schwarzschild-anti de Sitter line element with negative cosmological constant 
A = -3A 2 /4. 
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V. DISCUSSION 



We have highlighted here a curious aspect of McVittie's solution, namely that it displays 
an unstable aspect of the Schwarzschild event horizon and the black hole interior. Physical 



perturbations of the event horizon are governed by Price's theorem |TTJ and indicate its 



stability. However these perturbations do not extend to the situation considered here, where 
the black hole is immersed in an isotropic perfect fluid. The result is that the event horizon 
becomes singular, albeit weakly singular. An analogous situation obtains in the interior; 
the black hole interior cannot be embedded into an isotropic perfect fluid in the manner 
described above. It is therefore important to determine the following. Is this behaviour 
generic in any way ? That is, is McVittie's solution just one of a general class of solutions 
representing the physical embedding of the Schwarzschild space-time into a perfect fluid 
in which the event horizon becomes singular, or is this behaviour specific to McVittie's 
solution? Furthermore, one should ask if such models have any basis in reality; could such 
space-times be used to model matter accreting onto a black hole, for instance? With regard 
to the first question, removing the simplifying assumption of shear-free fluid flow, as would 
be necessary to discuss the more general situation, greatly increases the complexity of the 
field equations (see paper II for related comments). However the conditions described in 
section one above offer an indication of what boundary conditions should be used in solving 
Einstein's equations to address this question and these could lead to some simplifications. 

We have given an unphysical solution which applies to the embedding of the black hole 
interior into an isotropic medium. To obtain a physically (more) realistic solution, it would 
be helpful to find a geometric extension of the metric (|1.3|) through the weak space-like 
singular hypersurface R = 2m. The fact that the singularity is weak indicates that it may 
be possible to define, in some invariant or canonical way, a metric at R = 2m, using the 
full set of finite Jacobi fields at the surface. If it were possible to do the same for the first 
fundamental form at R = 2m, then it may be possible to set up a well posed initial value 
problem for the metric inside R < 2m. These ideas are very loosely stated, but the simple 
form of the singularity in this case may allow some progress to be made in this direction. We 
note that a good deal of progress has been made on a somewhat similar program which deals 
with space-times containing a different class of singularities which have a manageable degree 
of pathology, namely isotropic (or conformally compactifiable) singularities |T^]. For certain 
types of energy-momentum tensor, the conformal Einstein equations for such space-times 
have been shown to possess a well-posed Cauchy problem with initial data consisting of the 
conformal 3-metric at the singularity ||13| . 

We conclude by giving the latest version of the interpretation of McVittie's solution. The 
line element ( Jl.l] ) with u > m/2 gives the solution for the exterior region R > 2m of the 
Schwarzschild black hole embedded into a flat isotropic universe with scale-factor e^/ 2 ; the 
interior region cannot be embedded in this way. 
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